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ABSTRACT 

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It 
has been detected at radio frequencies and up to high-energy 7 rays (above 100 MeV). However, many models 
also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic 
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persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of 
the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hr 
between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of en- 
hanced 7-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An 
upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10~ 12 photons cm~ 2 s _1 
(95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection 
of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible 
VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The 
current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions. 
Subject headings: acceleration of particles — binaries: general — gamma rays: general — methods: observa- 
tional — X-rays: binaries 



1. INTRODUCTION 

Cygnus X- 3 is a bright and persi stent X-ray binary, discov- 
ered in 1966 dGiacconi et al.l ll967). but the high-energy pro- 
cesses occurring in the source are still poorly understood. It 
lies close to the Galactic plane at a distance between 3.4 and 
9.8 kpc, probably at 7 kpc, depending on differ ent distanc e 
estimates to the Cygnus OB2 association (Ling et al. 2009). 
The nature and the mass of the compact object are still a sub- 
ject for deb ate. Published result s suggest either a neutron star 
of 1.4 (IStark fc Saial 120031) or a black hole of less than 
10 Mq (Hanson et al.ll2000l). The identification of it s donor 
star as a Wolf-Rayet star (Ivan Kerkwiik et al.l [19921) classi- 
fies it as a high-mass X-ray binary. Nevertheless Cygnus X- 
3 shows a short orbital period of 4.8 hr, typical of the low- 
mass binaries, wh ich has been inferred fr om the modulation 
of bot h the X-rav dParsignault et al.ll 19721) and infrared emis- 
sions dBecklin et al.lll973l) . 

Despite the strong X-ray absorption in this system, which 
may be caused by the wind of the companion star, the X- 
ray spectrum has been intensively studied. The source shows 
two main spectral X-ray states resembling the canonical states 
of the black hole binaries: the hard state (HS) and the soft 
state (SS; Zdziarski & Gierlinski 2004; Hjalmarsdotter et al. 
2008). The HS is characterized by a weak soft thermal com- 
ponent and a strong non-thermal power-law emission peak- 
ing at hard X-ray e nergies, whereas th e SS, though showing 
a non-thermal tail (Szosteket al. 2008), is dominated by the 
optically thick thermal disk emission. In Cygnus X-3, how- 
ever, the HS displays a high-energy cutoff at rs 20 keV, signif- 
icantly lo wer than the ~ 100 keV va l ue found for black hole 
binar ies dHialmarsd otter et al.l I2004t [Zdziarski fc Gierlinskil 
12001) . 

Adding to its peculiarity, Cygnus X-3 is the brightest ra- 
dio source among the X-ray binaries, quite frequently ex- 
hibiting huge radio flares (as intense as few thousand times 
the quiescent emissi on level of ~ 20 mJy at 1 .5 GHz ; 
iBraes & Milevl 119721) first seen by iGregorv et all (119721) . 
During these outbursts, which occur mainly when the 
source is in the SS, and last several days, Cygnus X- 
3 reveals the presence of collimated relativistic jets (e.g.. 



Martf et al.ll200lUMioduszewski et al.H2001t fGeldzahl er et all 
19831; iMiller- Jones et al.ll2004l) . Thus. Cygnus X-3 has been 
classified as a microquasar. Based on arcsecond-scale ra- 
dio observations and their relation with soft X-ray emis- 
sion, six X-ray/radio states have been identified: quiescent 
(flux densities ~ 100 mJy), minor flaring ( < Uy), sup- 
pressed (< 100 mJy), que nched (< 30 mJy), major flaring 
(> 1 Jy), and post-flaring (fSzostek et alJl2008t iTudose et all 
120 lOt iKolionen et al.l 120 101) . In the transition from the X- 
ray hard/radio quiescent state to the SS, the radio emission is 



strongly suppressed, and if it reaches the que nched level, the 
source usually produces a major radio flare (Waltm an et al.l 

cm. 

Cygnus X-3 has also historically drawn a great deal 
of attention due to numerous claims of detection at TeV 
and PeV 7 rays, using early-days detectors in these en- 



1971 iDanaher et ail. fl98Tt 
)owthwaite 1983; Samorskv & Stamm 



ergy ranges jVladimirs kv et al 
Lam beFaLl 119821; ' 
1983[ ICawley et al. 

1986). However, a critical analysis of these observa 
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tions raised doubts about their validity (Chardi n fc Gerbier] 
1989), and in recent years the more sensitive instruments 
have not confirmed those claims for energie s above 500 
GeV (ISchilling et al.l l200Tt lAlbert et all l2008al) . Neverthe- 
less, microquasars are believed to produce a very high en- 
ergy (VHE) emission inside the jets: high density and mag- 
netic fields provided by the accretion disk and by the compan- 
ion star cr eate favorable conditions for effective production 
of 7 ravs dLevinson & BlandforHHH iRomero et all [2003: 
iBosch- Ramon et al. 2006). This radiation could have either 
an epi sodic nature due to the eject ion of strong radio-emitting 
blobs (Ato van fe Aharonianll999l) . generally occurring in the 
SS in the case of Cygnus X-3 or a (quasi) stationary character 
being o riginated in the persistent compact jet present during 
the HS (IBosch-Ramon et al.ll2006l) . 

Using data from the Energetic Gamma-Ray Experiment 
Telescope detector aboard t he Compton Gamma-Ray Ob- 
servatory, Mor i et al.l (119971) reported an average flux of 
(8.2 ± 0.9) x 10~~ 7 photons cm~ 2 s _1 at energies above 
100 MeV coming from the direction of Cygnus X-3. How- 
ever, no orbital modulation was detected in the signal, pre- 
cluding a solid association. The experimental situation in the 
high-energy region has been drastic ally changed by th e re- 
sults recentl y published by A GILE dTavani e t al. 200% and 
Fermi/LKT lAbdo et al.l 12009). AGILE detected the source 
in five different moments, for a couple of days each, four of 
them corresponding to the peak emissions shown in the de- 
tailed Fermi/LAil light curve (see below). The last detection, 
occurred in 2009 July, has not been published yet by the AG- 
ILE Collaboration (A. Bulgarelli et al., in preparation). On the 
other hand, FermifLPtf detected a clear signal from Cygnus 
X-3 above 100 MeV during two periods of enhanced activity 
lasting for several weeks and coinciding with the source being 
in the SS. The measured flux is variable and shows an orbital 
modulation, which confirms the origin of the signal from the 
microquasar. The AGILE flux level is comparable with the 
one of the FermilLKI flux peaks, as high as 2.0 x 10~ 6 pho- 



tons cm 



above 100 MeV. 



Observations of binary systems with imaging atmospheric 
Cherenkov telescopes (IACTs) in the VHE band have proven 
very fruitful in recent years, with the detection of the orbitally 
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modul ated 7-ray e mitters PSRB 1259-63 |Aharonian et al.l 
l2005ah . LS 5 039 dAharonian et alJ I2005U I2006D. and 



LS 1 +61 303 dAlbert et al.ll2006l l2008bl 120091: lAcciari et alJ 
120081 12009b [Anderhub et al. 2009j). However, these sys- 
tems may be different from Cygnus X-3 since the radio 
and high-energy radiation could be produced in all three 
sources by the interaction of the winds of a star and an orbit- 
ing pulsar dMaraschi fc Treves! 1 1 98 U iTavani fc Aronslll997t 
Dubusl l2006l) . Nevertheless. Cygnus X-3 GeV 7-ray modu- 
lation presents some common features to those observed in 
LS5039 and LS I +61 303. On the other hand, although 
observations of ot her well-established microquasars, s uch as 
GRS 1915+105 (lAcero et all 120091 iSaito et al.l l2009h . did 
not reveal any signal, there is an interesting possibi lity of 
VHE emission from Cygnus X-l dAlbert et al.l l2007ah . still 
to be confirmed by an independent detection. The origin of 
this possible emission has been speculated to be linked to a 
maxi mum of the X-ray super-orbital modulation of the sys- 
tem dRicol 12008). Moreover, there is also a claim by AGILE 
about a Cygnus X-l detection during a short flare (Sabatini et 
al. 2010, ATel 2512), which, however, has not been corrobo- 
rated by Fermi/LAT yefl. Such experimental results endorsed 
by the theoretical predictions, have encouraged deeper obser- 
vations of the black hole binaries in the VHE band. 

This paper reports observations of Cygnus X-3 performed 
with the Major Atmospheric Gamma Imaging Cherenkov 
(MAGIC) telescope between 2006 and 2009. All these ob- 
servations were carried out with the first stand-alone MAGIC 
telescope, MAGIC phase I. Cygnus X-3 observations were 
planned to cover different X-ray/radio states, including those 
that showed a strong high-energy 7-ray flux. This allowed us 
to search for VHE emission, above 250 GeV, from Cygnus X- 
3 in the different cases for which this radiation is predicted in 
theoretical scenarios. In addition, specific analyses were per- 
formed to look for predicted features, such as periodic and 
variable emission. Section |2] describes the performance of 
MAGIC and the observational strategy. The analysis chain 
is explained in Section [3] The results obtained by MAGIC in 
a multi-wavelength context are reported in Section|4] A brief 
discussion is given in Section|5] 

2. OBSERVATIONS 

MAGIC (phase I) is an IACT located at the Canary Island 
of La Palma (Spain), at 28°.8 N-, 17°.8 W, 2200 m above sea 
level. It is a 17 m diameter IACT with an energy threshold of 
60 GeV (with the standard trigger). The collected Cherenkov 
light is focused on a multi-pixel camera composed of 576 pho- 
tomultiplier tubes (PMTs). 

The performance of the MAGIC telescope changed over the 
years. The largest improvement followed the upgrade of the 
300 MHz readout system in 2007 February. The new mul- 
tiplexed 2 GHz flash analog-digital converters improved the 
time resolution of the recorded shower images a nd reduced 
the c ontamination of the night sky background dAliu et alJ 
2009). Accordingly, the telescope integral flux sensitivity im- 
proved from ps 2.5% to « 1.6% of the Crab Nebula flux 
for energies greater than 270 GeV in 50 hr of observations. 
At these energies, the angular and energy resolutions are 
ps 0°.l and ps 20%, respectively. MAGIC can provide 7-ray 
point- like source localiza tion in the sky with a precision of 
pa 2' dAlbert et al.ll20 08c). It is able to observe under mod- 
erate moonlight or twilight conditions (Alb ert et al.l l"2007b: 



Preliminary results can be found at http://fermiskyMogspot.com/2010/03Aat-limiamL-cyg-x-l-during-reported.html 



Britzger et al. 2009), which, causing an increase of the night 
sky background, can be monitored through the PMT anode 
direct currents (DCs). 

MAGIC pointed toward Cygnus X-3 for a total of 69.2 hr 
between 2006 March and 2009 August. Since the source is 
expected to be variable, the observations were triggered by 
its state at other wavelengths. In 2006, observations were 
prompted by flares at radio frequencies. The MAGIC collab- 
oration received two alerts from the RATAN-600 telescope 
on 2006 March 10 and on 2006 July 26 (S. Trushkin, private 
communication 2006). Both radio flares had an X-ray coun- 
terpart: the source was in the SS according to RXTEI ASM 
and Swift/BAT measurements. In 2007, the observations 
were triggered using public RXTE/ASM (1.5-12 keV) and 
Swift/BAT (15-50 keV) data, as follows: (1) Swift/BAT daily 
flux larger than 0.05 counts cm~ 2 s _1 and (2) ratio between 
RXTE/ASM and Swift/BAT count rates lower than 200. These 
criteria guaranteed the source to be in the HS. During 2008 
and 2009, MAGIC observed Cygnus X-3 following high- 
energy 7-ray alerts issued by the AGILE satellite. The first 
of these alerts arrived in 2008 April 18, the second one in 
2009 July 18 (M. Tavani, private communication 2008; 2009). 
These last two campaigns occurred when the source was in the 
SS. 

At La Palma, Cygnus X-3 culminates at a zenith angle 
of 12°. The observations were carried out at zenith angles 
between 12° and 50°. They were performed partially in 
a on -off mode and in a false-source tracking mode (wob- 
ble) dFomin et al Jl 19941) . the latter pointing to two directions 
at 24' distance and opposite sides of the source. The entire 
Cygnus X-3 data set recorded by MAGIC amounts to 69.2 
hr, out of which 12.5 hjQ were rejected from further analysis 
because of their anomalous event rates, leading to a total of 
56.7 hr of useful data distributed in 39 nights of observation. 
The details of the observations are quoted in Table Q] 



3. DATA ANALYSIS 

Data analysis was carried out using the standard MAGIC 
calibration and analysis software: o nce the PMT signal pulses 
are calibrated (Albert et al. 2008d), pixels containing no use- 
ful information for the shower image reconstruction are dis- 
carded. This is done by an image cleaning procedure which 
takes into account the amplitu de and the timing information 
of the calibrated signals dAliu e t al. 2009). The constraints on 
the signal amplitude are increased in the case of moonlight 
conditions, when the pixel DCs become larger than 2.5 /LtA 
(see Table Q]). This higher image cleaning allows us to use 
the standard Monte Carlo (MC) during the analysis of these 
data, without the need of any c orrection, as long a s the aver- 
age pixel DCs are below 4 fj,A dBritzger et alJl2009h . 

Afterward, the surviving pixels are used to compute the 
Hillas event image parameters (Hillai 119851) . In addi- 
tion, Hillas and timing variables are combined into a sin- 
gle 7/hadron discriminator (hadronness) and an energy esti- 
mator by means of the Random Forest classification and re- 
gression algorithm, which takes into account the c orrelation 
between the different Hi llas and timing variables dBreimanl 
1200 HI Albert et alj|2008el) . These algorithms are trained with 
a sample of MC-simulated 7-ray events and a sample of back- 

2 9.6 hr were rejected due to high-altitude Saharan dust (calima), which 
affects the atmosphere above Canary islands and is more intense during sum- 
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TABLE 1 
Cygnus X-3 Observation log 3 
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a From left to right: observational cycle, date of the beginning of the 
observations, also in MJD, effective time after quality cuts, zenith angle 
range, the anode PMT DCs, and the spectral state. 

b The anode PMT DCs show the level of the moonlight conditions. 

c Spectral state was derived by using SwiftfBAI and RXTE/ASM data. 



with respect to the camera center. 

The Cygnus X-3 data set extends over five different, one- 
year long cycles of observation, which are characterized by 
varying performances of the telescope. In addition, cycle 
I data were recorded in the on-off mode, all the others in 
the wobble mode, cycle IV data were taken under moderate 
moonlight conditions. 

Each cycle of data was analyzed separately using the ap- 
propriate image cleaning procedure, and a matching sample 
of MC-simulated 7-ray events. The analysis selection cuts, 
on hadronness and alpha, for each cycle of data were obtained 
optimizing the sensitivity in a Crab Nebula sample and requir- 
ing at least 70% 7-ray selection efficiency. The selected Crab 
Nebula sample was recorded during the same cycle at similar 
zenith angles and in the same data taking mode. 

All the analyses were then combined in order to calculate 
the integral flux upper limits (ULs) for energies greater than 
250 GeV, which is the energy threshold of each cycle of data. 
The obtained ULs at the 95% co nfidence level (CL) were 
computed after Rol ke et all (|2005) method assuming a Gaus- 
sian backg round and 30% of systematic uncertainties in the 
flux level (Alb ert et alj|2008c|). The sp ectrum was assumed 
to be Crab like ( Aharoni an et alj |2004). with a photon index 
of 2.6. However, a 30% change in the photon index yields a 
variation of less than 1% in the ULs. 

4. RESULTS 

A search for VHE 7-ray emission from Cygnus X-3 was 
performed separately for each cycle of observations. None of 
them yielded a significant excess. The combination of all the 
data samples yields a 95% CL upper limit to an integral flux 
of 2.2 x 10~ 12 photons cm~ 2 s _1 for energies greater than 
250 GeV. It corresponds to 1.3% of the Crab Nebula flux at 
these energies. The differential flux ULs are shown in Table|2] 
and in Figure [TJ 

Given that Cygnus X-3 flux is variable at other energy 
bands on time scales of days, 7-ray signals were searched for 
also on a daily basis. No significant excess events were found 
in any observation night. The integral flux ULs for energies 
above 250 GeV are shown in Tableland in the top panel of 
Figure |2 



ground events extracted from real data. The used timing vari- 
ables allow an improvement of the analysis sensitivity of a 
factor of 1 .4 ( 1 .2) in data re corded with the n ew ultra-fast (old 
300 MHz) readout system dAliu et alj|2009l) . 

Images with a total charge below 200 photo-electrons were 
discarded from further analysis in order to homogenize the 
different data and MC samples and achieve a common stable 
energy threshold of all the analyses. The number of 7-ray 
candidates in the direction of th e source is e stimated by using 
the distribution of alpha angle (Hillas 1985), which is related 
to the shower orientation. The evaluation of the background 
depends on the data taking mode. For data taken in the on- 
off mode in which the signal data sample is called "on" data, 
the background is estimated using a different sample, called 
"off data. The latter has a similar amount of data recorded 
during the same epoch as the "on" sample at similar zenith an- 
gles and atmospheric conditions with no known 7-ray source 
in the field of view. For the observations carried out in the 
wobble mode, the background is extracted from three circular 
control regions, located symmetrically to the source position 




FIG. 1. — Differential flux ULs at 95% CL for the VHE time-integrated 
emission. The slope of the arrows indicates the assumed power-law spectrum 
with a photon index of 2.6. 
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FIG. 2.— From top to bottom: daily MAGIC VHE integral flux ULs for E > 250 GeV, high-energy 7-ray (AGILE and Fermi/LAT), hard X-ray (Swift/BAT), 
soft X-ray (RXTEI ASM), and radio fluxes measured for Cygnus X-3 as a function of time (from 2006 January 1 until 2009 December 15). The grey bands show 
the periods corresponding to the MAGIC observations. The OVRO and AMI 15 GHz data generally agree well, except for the offset apparent during steady 
periods which is due to unrelated extended emission resolved out by AMI. 



In Figure [2] MAGIC results are presented in a multi- 
wavelength context where the above-mentioned flux vari- 
ability at different energy bands is rather clear. In par- 
ticular, it displays daily MAGIC VHE integral ULs, high- 
energy 7-ray (AGILE and FermifLAT (0. 1-30 GeV)), hard X- 
ray (Swift/BAT (15-50 keV)), soft X-ray (RXTEI ASM (3-5 
keV)), and radio measured fluxes from 2006 January 1 (MJD 
53736) until 2009 December 15 (MJD 55180). The soft X- 
ray energy band of RXTEI ASM is between 1.5 and 12 keV, 
out of which only the 3-5 keV band was used, as it yields to 



the cleanest radio/X-ray correlation dSzostek et al.l2008h . The 
radio measurements, displayed in logarithmic scale, were pro- 
vided by the RATAN-600 telescope at 2.15, 4.8, and 1 1 .2 GHz 
and by the Ryle telescope (RT), the Arcminute Microkelvin 
Imager (AMI) telescope and the Owens Valley Radio Obser- 
vatory (OVRO) 40 meter telescope at 15 GHz. 

The soft and hard X-ray fluxes shown in Figure [2] allow 
us to derive the X-ray spectral state of Cygnus X-3 during 
MAGIC observations. Cycle III data are the only ones taken 
when the source was in the HS, as requested by the obser- 
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FIG. 3. — Zoom of Figure [2] around the cycle V campaign between 2009 May 29 and August 17. The open black points in the VHE MAGIC panel show the 
non-significant measured integral fluxes with their statistical error bars (whereas the ULs take into account also the systematic errors). 



vational trigger. All the other observations were carried out 
when the source was in the SS. The first two MAGIC obser- 
vational campaigns (cycles I and II) were triggered by flares 
at radio frequencies, which are expected when the source is in 
the SS. Unfortunately, the conditions on the trigger together 
with MAGIC constraints on observation time did not allow a 
simultaneous coverage of the flaring states: MAGIC started 
pointing toward Cygnus X-3 twelve and two days after the 
strong radio emission, respectively. On the other hand, in the 
last observational cycles (IV and V), MAGIC observed the 
source during the SS, following two AGILE alerts on a high- 
flux activity in the high-energy band: the GeV emission oc- 



curs only in this X-ray spectral state. During the second alert, 
in 2009 July (MJD 55030), MAGIC recorded some data si- 
multaneous with a GeV peak emission and did not detect any 
VHE sign al. T his important result will be discussed in detail 
in Sectionl4~T1 

Microquasars are expected to produce VHE emission inside 
the jets either when they are compact and persistent, mainly 
in the HS, or in the presence of strong radio-emitting-blobs, 
which most likely happen during an SS. Therefore, the re- 
sults of the MAGIC observations were divided according to 
the X-ray sp ectra l state of the source, as described in Sec- 
tions [4j2] and [43] for the SS and HS, respectively. Besides, 
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in the scenario where the VHE radiation is emitted during 
powerful synchrotron radio ejections, th ere might be a corre- 
lation between these two w avele ngths dAtovan fc Aharonianl 
1999). Therefore, in Section l4~4l the X-ray/radio states during 
the MAGIC observations are quoted. 

4. 1 . Results during High-energy "/-ray Emission 



lAbdo et a l. (2009) presented the FermifLPS observations of 
Cygnus X-3 between 2008 August 4 and 2009 September 2. 
They detected a strong signal above 100 MeV, with an over- 
all significance of more than 29 standard deviations, which 
is mainly dominated by two active periods: MJD 54750- 
54820 and 54990-55045. These active periods coincide with 
the X-ray SS of the source, indicating that Cygnus X-3 emits 
high-energy 7-rays during this spectral state, with an average 
flux of 1.2 x 10 -6 photons cm~ 2 s _1 and a photon index of 
2.70 ± 0.05 s tat i 0.20 sys (under the assumption of a power- 
law spectrum). They also estimated that the peak flux can be 
as high as 2 x 10~ 6 photons cm~ 2 s _1 without providing any 
photon index. 

The five AGILE detections between 100 MeV and 3 GeV 
(Tavani et al. 2009; A. Bulgarelli et al., in preparation) co- 
incide with the strongest high-energy outbursts, which can be 
seen overlapping them with the FermilhPC[ light curve. The 
AGILE integral flux averaged over the first four detections is 
estimated to be 1.9 x 10~ 6 photons cm~ 2 s _1 with a corre- 
sponding photon index of 1.8 ± 0.2. 

In cycle V, MAGIC pointed at Cygnus X-3 during the 
second period of high-energy enhanced activity detected by 
fen«//LAT, as shown in Figure [3] In particular, MAGIC car- 
ried out observations simultaneous with a GeV emission peak 
on 2009 July 21 and 22 (MJD 55033-55034), and it did not 
detect any VHE emission. The corresponding MAGIC inte- 
gral flux ULs above 250 GeV are lower than 6% of the Crab 
Nebula flux. 

Figure |4] shows the spectral energy distribution (SED) of 
Cygnus X-3 between 100 MeV and 5 TeV including MAGIC 
95% CL upper limits at VHE, and the power-law spec- 
trum in the high-energy range reported by both AGILE and 
Fermi/LAiT. The spectra take into account the error on the 
photon index and the one on the integral flux, which is 30% 
and 40% for FermifLKT and AGILE, respectively. The SED 
for the average SS was obtained considering the average 
Fermi/LAn: flux and MAGIC results of the SS data set. On 
the other hand, to obtain the SED during a high-energy emis- 
sion peak, the MAGIC measurements simultaneous with the 
GeV emission peak and both Fermi/LAT and AGILE spectral 
power-law fits were used. The FermifLPZI photon index dur- 
ing the peak was assumed to be the same as the one for the 
SS average flux. Being the latter dominated by the bright- 
est phases of the 7-ray outburst, it can also be considered 
representative of the flaring activity. Both MAGIC ULs are 
consistent with the extrapolation of the fermZ/LAT spectra up 
to VHE, but not with the extrapolation of the harder AGILE 
spectrum, which would suggest a cutoff between some tens 
and 250 GeV. 

4.2. Results During the Soft State 

The MAGIC telescope pointed at Cygnus X-3 when it was 
in the SS during the observational cycles I, II, IV and V, corre- 
sponding to a total of 30.8 hr. For these observations, soft X- 
ray measurements in the 3-5 keV band are always above th e 
transitional level set at 3 counts s _1 by Szos tek et alj ([2008). 
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After having analyzed each data cycle separately, the data 
sets were combined in order to obtain a global UL to the in- 
tegral flux for the SS. The UL at 95% CL for energies greater 
than 250 GeV is 4.1 x 10~ 12 photons cm" 2 s _ \ i.e., ~ 2.5% 
of the Crab Nebula flux. The differential flux ULs for this 
spectral state are shown in Table [4] and in the left panel of 
Figure [5] 

Due to the highly anisotropic radiation from the compan- 
ion star, the predicted 7-ray emission above 250 GeV woul d 
be modulated according to the orbital phase (Bednarek 1997). 
Absorption might play an important role in making the VHE 
orbital modulation difficu lt to be detecte d by the current sen- 
sitivity instrumentation (Bednarek 2010). The orbital modu- 
lation of the GeV 7-ra y emission was de tected only when the 
source was in the SS (lAbdo et alJl2009l) . MAGIC searched 
for such modulation in this spectral state. A phase-folded 
a nalysis was perform ed assuming the parabolic ephemeris 
in lSingh et al.l ((2002). The results are shown in the left panel 
of Figure [6] and in Table [5] No evidence of VHE 7-ray signal 
was found in any phase bin. The obtained integral flux ULs 
are smaller than 10% of the Crab Nebula flux for all of them. 

4.3. Results During the Hard State 

The 25.9 hr of MAGIC cycle III data sample were ob- 
tained when Cygnus X-3 was in the HS. Svw/f/BAT count 
rates during this cycle are rather high, greater than 0.05 
counts cm~ 2 s _1 , whereas the soft X-ray fluxes in the 3-5 
keV band are below 2 counts s _1 . 

For this spectral state, the integral flux UL for energies 
greater than 250 GeV is 1.1% of the Crab Nebula flux 
(1.8 x 10~ 12 photons cm~ 2 s _1 ). The differential flux ULs 
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FIG. 5. — Ninety-five percent differential flux ULs for the SS (left panel) and HS (right panel) observations. The slope of the arrows indicates the assumed 
power-law spectrum with a photon index of 2.6. 
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are quoted in Tableland plotted in the right panel of Figure[5] 
The performed phase-folded analysis for this spectral state did 
not yield any significant detection. The computed ULs to the 
integral flux are, on average, at the level of 5% of the Crab 
Nebula flux (see the right panel of Figure|6]and Table|5]l. 

4.4. Results During X-ray/Radio Rtates 

MAGIC observed Cygnus X- 3 in both X-ray main spec- 
tral states (see Sections |4.2| and |4~3l . However, the state of 
the s ource can be further c haracterized by simultaneous radio 
flux. iSzostek et al.l (|2008) identified six different X-ray/radio 
states studying simultaneous observations of the Green Bank 
Interferometer at 8.3 GHz and RXTE/ASM in the energy 
range 3-5 keV. The relation between these two energy bands 
is shown in the so-called saxophone plot (see Figure|7]i. It was 



noted that the use of other radio frequencies yields similar re- 
sults. This gives us the confidence that a direct comparison 
between their and our results using 15 GHz (RT and AMI) 
and 11.2 GHz (RATAN-600) is reasonable. The OVRO and 
AMI 15 GHz data generally agree well, but for an ^0.12 Jy 
offset apparent during steady periods, probably due to unre- 
lated extended emission resolved out by AMI. Thus, only the 
AMI 15 GHz were used in this analysis, although our conclu- 
sions are not substantially affected by this choice. 

Figure|7]shows the soft X-ray (3-5 keV) RXTE/ASM count 
rates versus radio flux densities corresponding either to the 
nights of MAGIC observations for the five observational cy- 
cles or to the AGILE flux peaks (only the last four AGILE 
detections were considered since no simultaneous radio data 
were available for the first one). Unfortunately, no radio 
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data simultaneous with cycle III and cycle IV MAGIC ob- 
servations were available. Nevertheless, in the case of cy- 
cle III data, this does not prevent us from pointing out that 
Cygnus X-3 was in a quiescent state, just by using the soft 
X-ray measurements. For cycle IV data, quasi simultane- 
ous radio observations (one day before MAGIC observations) 
were used. This latter choice does not affect our qualita- 
tive result since the source had already entered a post-flaring 
state. As shown in Figure [H] in 2008 April, AGILE detected 
Cygnus X-3 (MJD 54572-545 73) one day before a major ra- 
dio flare dTrushkin et alj20 08) lasting a few days, but MAGIC 
started pointing at the microquasar 10 days after the peak ra- 
dio emission. 

All the high-energy flux peaks were detected in the right 
branch of the " saxophone" , befor e, after, or during a flaring 
state. Although lAbdo et aD ([2009) quoted a time lag between 
the radio and the 7-ray peaks of 5 ±7 days, the correlation 
between the two energy bands is not yet clear. On the other 
hand, MAGIC observed Cygnus X-3 in its SS some days after 
the radio flare occurred (see Figures [8] 0, although for the 
first nights of cycle V observations, the radio flux densities 
are rather high and oscillating between two small flares (see 
Figure [3]). 

5. DISCUSSION 

MAGIC observations of Cygnus X-3 cover all X-ray spec- 
tral states of the source in which VHE emission is thought to 
be likely produced either from a persistent jet in the HS or 
during powerful ejections in the SS. However, no significant 
excess events were found in any of the inspected samples. 

VH E 7 rays have been predi ct ed from micr o quasar 
jets ( lAtovan fe Aharonianl Il999t iRomero et all 120031: 



iBosch-Ramon et ai1 f2006). A robust prediction of modeling 
is that photon-photon absorption cannot be ne glected if 
the 7 rays are produ ced close to a massive star dBednarekl 
119971: lOrellana et ail 120071) . In particular for Cygnus X-3 
the presence of a Wolf-Rayet companion, with temperature 
T* ~ 10 5 K and radius i?„ 10 11 cm, leads to an optical 
depth t > 1 for VHE 7 rays for an em i tter lo cated at several 
orbital radii from the star (Bednarek 2010). Even under 
very efficient electromagnetic cascading, i.e., a radiation to 
magnetic energy density ratio 8iru*/B 2 1, the expected 
VHE fluxes are below the sensitivity of the present instru- 
ments (Be dnarekll2010l) . Therefore, in order to detect VHE 
photons, the emitter should be located far from the binary 
system. 

Fermi/LAI detected Cygnus X-3 when it was in the SS and 
found orbital modulation for the radiation above 100 MeV 
with a photon index of 2.7 for the periods of enhanced ac- 
tivity (Ab do et al.ll2009l) . For the epochs outside these high- 
activity pe riods, the GeV fl ux decreases significantly (Fig- 
ure |2] in lAbdo et ail (120091) ) and no modulation is found. 
The GeV orbital light curve of Cygnus X-3 in the high- 
energy active periods can be explained in the context of 
anisot ropic inverse Compton scattering with the stellar pho- 
tons (lAbdo et al.l|2009b IBosch-Ramon fc Khangulyanl 120091: 
iDubus et al.ll2010l) . which is also energetically more efficient 
than hadronic mechanisms such as pp interactions or pho- 
tomeson production. Since very bright X-ray emission is pro- 
duced in the inner accretion disk or at the base of the jets in 
Cygnus X-3, the GeV radiation would be absorbed unless it 
is originated beyond ~ 10 10 cm above the compact object. 
This implies that the GeV radiation is produced in the jet of 
Cygnus X-3 rather than in the inner accretion disk/corona re- 
gion. On the other hand, the GeV emitter cannot be too high 
in the jet, since otherwise there would not be strong orbital 
modulation (see also Abdo et al. 2009). Therefore, the ob- 
served GeV and the predicted detectable VHE emission can- 
not be explained by one particular population because the for- 
mer should be produced in/close to the system and the latter 
farther from it. The location of a hypothetical VHE emitter 
could coincide with the innermost region of the radio emitting 
jet, which, to avoid synchrotron self-absorption, should start 
relatively far from the binary system. This is consistent with 
the fact that MAGIC did not detect Cygnus X-3 during the 
high-activity GeV period, even though the flux ULs are close 
to a power-law extrapolation of the FermilhP^L spectrum to 
energies greater than 100 GeV and well below an extrapola- 
tion of the AGILE spectrum. 

During the periods of GeV high-activit y of the source, as by 
the sy nchrotron self Compton scenario ( lAtovan fc Aharonianl 
1999), a detectable TeV signal could arise during the first 
hours of a radio outburst. Unfortunately, MAGIC has never 
observed the source during this phase of the flare, but always 
some days before or after the maximum radio flux. This ra- 
diation would not be strongly modulated, due to its origin far 
from the system. The two times more sensitive two telescopes 
arrangement, MAGIC phase II, may indeed detect Cygnus X- 
3 if it observes the source for longer time at the very maximum 
of a GeV flare, w hich might be even earlier than the onset of 
the radio outburst (I Abdo et al.ll2009l) . 

In the HS, the VHE emission is expected to be produced 
inside the compact and persistent jets, who se total luminos - 
ity is estimated to be at least 10 37 erg s -1 dMartf et al.ll 2005). 
The MAGIC VHE 7-ray UL set at 1.1 x 10" 12 erg cm" 2 s" 1 
is equivalent to a VHE luminosity of ~ 7 x 10 erg s _1 
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at 7 kpc. Thus, the maximum conversion efficiency of 
the jet power into VHE 7 rays is 0.07% which is sim- 
ilar to that of Cygnus X-l for the UL on the VHE 
steady emission, bu t 1 order of magnitude greater than tha t 
of GRS 1915+105 dAlbert et all l2007at lAcero et all l2009l) . 
These ULs are in good agreement with the theoretical expec- 
tations which generally predict a VHE steady luminosity of ~ 
10 32 erg s -1 . Persistent galactic jets do not seem to be good 
candidate sources to be detected at VHEs by the current sen- 
sitivity instrumentation. Only 10 times more sensitive future 
instruments, such as Cherenkov Telescope Array, may have a 
chance to detect such VHE emission. This would provide a 
new handle on the emission mechanisms of compact jets. 
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MAGIC constraints on 7-ray emission from Cygnus X-3 1 1 



TABLE 2 

Differential Flux ULs for the VHE Time-integrated Emission at 95% CL. 



Energy Range 


N on EvtS." 


N bg Evts. b 


Excess Evts. 


Norm.Fact. c 


Signif. d 


UL Evts. c 


Flux UL 


(GeV) 










(T) 




(TeV- 1 cm- 2 s- 1 ) 


199-315 


4416 


4384.5 ± 39.0 


31.5 ± 77.0 


0.34 


0.4 


237.8 


2.6E-11 


315-500 


2057 


1980.6 ±28.6 


76.4 ± 53.6 


0.36 


1.5 


264.2 


1.2E-11 


500-792 


769 


800.8 ±21.1 


-31.8 ±34.9 


0.39 


-0.9 


51.3 


LIE- 12 


792-125 


289 


299.9 ± 12.9 


-10.9 ± 21.4 


0.38 


-0.3 


39.2 


5.1E-13 


1256-1991 


102 


98.3 ± 6.9 


3.7 ± 12.2 


0.37 


0.5 


36.2 


3.0E-13 


1991-3155 


38 


32.3 ± 3.5 


5.7 ± 7.1 


0.35 


0.7 


27.4 


1.3E-13 



a Number of signal events. 



b Number of normalized background events. 
c Normalization factor. 
d Significance. 

c Maximum number of excess events computed by using Rolke's method. 



TABLE 3 

Integral Flux ULs for Energies Above 250 GeV Calculated on a Daily Basis at 95% CL. 



Date 


Time 


Non EvtS. 


N bg Evts. 


Excess Evts. 


Norm. Fact. 


Signif. 


UL Evts. 


Flux UL 




(MJD) 


(h) 










(t) 




(cm" 2 s -1 ) 


(% C.U.) 


53817 


0.46 


42 


41.4 ± 1.8 


0.6 ± 6.7 


0.07 


0.1 


18.7 


1.2E-11 


7.5 


53818 


0.26 


43 


27.0 ± 0.7 


16.0 ± 6.6 


0.05 


2.8 


46.7 


5.7E-11 


34.2 


53820 


0.54 


85 


83.8 ±4.6 


1.2 ± 10.2 


0.15 


0.1 


27.4 


1.5E-11 


9.3 


53822 


0.70 


132 


111.4 ±3.5 


20.6 ± 12.0 


0.19 


1.7 


67.9 


2.9E-11 


17.4 


53824 


0.80 


84 


84.8 ± 2.4 


-0.8 ± 9.5 


0.15 


-0.1 


24.3 


9.1E-12 


5.5 


53825 


0.80 


83 


64.3 ±2.1 


18.7 ±9.3 


0.15 


2.1 


58.0 


2.2E-11 


13.1 


53826 


1.00 


37 


34.3 ±0.6 


2.7 ± 6.1 


0.06 


0.4 


20.9 


6.1E-12 


3.7 


53827 


0.92 


79 


76.5 ± 1.9 


2.5 ± 9.1 


0.13 


0.3 


28.4 


9.1E-12 


5.5 


53828 


1.05 


42 


49.3 ± 1.1 


-7.3 ± 6.6 


0.09 


-1.0 


10.9 


3.1E-12 


1.8 


53943 


3.10 


257 


274.3 ± 9.5 


-17.3± 18.6 


0.33 


-0.9 


27.9 


4.1E-12 


2.5 


53944 


2.54 


232 


235.7 ± 8.8 


-3.7 ± 17.6 


0.33 


-0.2 


38.6 


7.1E-12 


4.3 


53945 


1.73 


122 


125.7 ± 6.4 


-3.7 ± 12.7 


0.33 


-0.3 


27.1 


8.4E-12 


5.1 


53946 


1.37 


127 


108.0 ± 5.9 


19.0± 12.7 


0.33 


1.5 


65.0 


2.5E-11 


15.3 


53948 


0.90 


82 


65.7 ± 4.6 


16.3 ± 10.2 


0.33 


1.7 


54.3 


3.2E-11 


19.4 


53949 


1.79 


138 


143.3 ± 6.9 


5.3 ± 13.6 


0.33 


0.4 


42.1 


1.3E-11 


7.8 


54286 


2.16 


299 


314.0 ± 10.2 


-15.0 ±20.1 


0.33 


-0.74 


62.1 


1.1E-11 


6.8 


54287 


4.53 


547 


547.0 ± 13.4 


-0.7 ± 27.0 


0.33 


-0.03 


150.3 


1.3E-11 


8.1 


54288 


1.07 


176 


208.0 ± 8.3 


-32.0 ± 15.6 


0.33 


-1.99 


20.3 


7.0E-12 


4.3 


54289 


0.38 


86 


84.3 ± 5.3 


1.7 ± 10.7 


0.33 


0.16 


41.3 


3.9E-11 


23.5 


54295 


1.82 


189 


176.3 ± 7.6 


12.7 ± 15.7 


0.33 


0.82 


76.4 


2.1E-11 


13.0 


54296 


1.93 


219 


209.0 ± 8.3 


10.0 ± 17.0 


0.33 


0.59 


74.9 


2.0E-11 


11.8 


54297 


1.93 


151 


138.7 ± 6.7 


12.3 ± 14.0 


0.33 


0.89 


65.9 


1.7E-11 


10.4 


54305 


1.75 


179 


167.7 ± 7.4 


11.3 ± 15.3 


0.33 


0.75 


111.6 


3.1E-11 


19.0 


54307 


0.98 


115 


117.0 ±6.2 


-2.0 ± 12.4 


0.33 


-0.16 


52.9 


2.4E-11 


14.5 


54308 


0.5 


71 


59.3 ±4.4 


11.7 ±9.5 


0.33 


1.27 


100.1 


8.1E-11 


48.8 


54316 


0.73 


85 


88.7 ± 5.4 


-3.7 ± 10.7 


0.33 


-0.34 


33.0 


1.8E-11 


11.0 


54318 


2.72 


262 


277.3 ± 9.6 


-15.3 ± 18.8 


0.33 


-0.81 


39.2 


6.7E-12 


4.1 


54319 


1.93 


184 


184.0 ±7.8 


0.0 ± 15.6 


0.33 


0.00 


61.8 


1.5E-11 


8.8 


54320 


1.58 


146 


152.3 ±7.1 


-6.3 ± 14.0 


0.33 


-0.95 


38.8 


1.1E-11 


6.8 


54346 


1.86 


182 


200.0 ±8.1 


-18.0 ± 15.7 


0.33 


-1.12 


45.3 


1.0E-11 


6.3 


54584 


0.31 


17 


23.0 ± 2.7 


-5.0 ±5.0 


0.33 


-1.0 


8.0 


1.0E-11 


6.2 


54585 


1.07 


89 


75.3 ±5.0 


13.7 ± 10.7 


0.33 


1.3 


50.2 


2.0E-11 


11.9 


54586 


1.33 


60 


66.7 ± 4.7 


-6.7 ±9.1 


0.33 


-0.7 


15.5 


5.1E-11 


3.1 


55031 


3.48 


186 


183.3 ±7.8 


2.7 ± 15.7 


0.33 


0.2 


43.2 


5.7E-12 


3.4 


55033 


2.21 


71 


63.7 ±4.6 


7.3 ± 9.6 


0.33 


0.8 


36.2 


7.2E-12 


4.4 


55034 


1.63 


50 


43.3 ±3.8 


6.7 ± 8.0 


0.33 


0.9 


31.4 


9.7E-12 


5.8 


55044 


1.81 


88 


80.0 ±5.1 


8.0 ± 10.7 


0.33 


0.8 


40.0 


9.4E-12 


5.7 


55045 


0.88 


69 


69.0 ± 4.8 


0.0 ± 9.6 


0.33 


0.0 


24.2 


1.1E-11 


6.6 



Note. — 



Refer to Table[2]for the meaning of the columns. 



12 



Aleksic et al. 



TABLE 4 

Differential Flux ULs for the SS and HS Observations. 



Spectral 


Energy Range 


N on Evts. 


N bo Evts. 


Excess Evts. 


Norm. Fact. 


Signif. 


UL Evts. 


Flux UL 


State 


(GeV) 








(a) 




(TeV _1 cm" 2 s _1 ) 




199-315 


1709 


1677.54 ± 24.6 


31.5 ±48.1 


0.36 


0.7 


169.2 


3.5E-11 




315-500 


926 


858.3 ± 18.5 


67.7± 35.6 


0.40 


1.9 


212.3 


1.7E-11 


HS 


500-792 


324 


357.5 ± 12.9 


-33.5 ± 22.2 


0.47 


-1.1 


30.3 


1.2E-12 




792-1256 


125 


124.9 ± 7.7 


0.1 ± 13.6 


0.48 


0.4 


38.0 


8.4E-13 




1256-1991 


41 


33.7 ± 3.9 


7.3 ±7.5 


0.46 


1.4 


32.7 


4.5E-13 




1991-3155 


14 


8.7 ± 1.8 


5.3 ±4.2 


0.40 


1.3 


20.3 


1.7E-13 




199-315 


2707 


2707.0 ± 29.9 


0.0 ± 60.0 


0.33 


0.0 


146.2 


3.4E-11 




315-500 


1131 


1122.3 ± 19.2 


8.7 ± 38.7 


0.33 


0.2 


108.5 


1.1E-11 


LH 


500-792 


445 


443.3 ± 12.1 


1.7 ± 24.3 


0.33 


0.1 


62.5 


3.3E-12 




792-1256 


164 


175.0 ± 7.5 


-11.0 ± 14.9 


0.33 


-0.7 


24.7 


7.7E-13 




1256-1991 


61 


64.7 ± 4.6 


-3.7±9.1 


0.33 


-0.4 


18.4 


3.7E-13 




1991-3155 


24 


23.7 ± 2.7 


0.3 ±5.6 


0.33 


0.1 


15.2 


1.7E-13 



NOTE. — Refer to Table[2]for the meaning of the columns. 



TABLE 5 

Integral Flux ULs for Energies Above 250 GeV for the Phase-folded Analyses of the Observations in the SS and the HS. 



Spectral 


Phase 


Time 


N on Evts. 


N bg Evts. 


Excess Evts. 


Norm. Fact. 


Signif. 


UL Evts. 


Flux UL 


State 




(h) 














(cm 2 s 1 ) 


(% C.U.) 




0.0-0.1 


1.34 


64 


68.0 ±4.7 


-4.3 ± 9.3 


0.33 


-0.5 


18.2 


6.7E-12 


4.0 




0.1-0.2 


1.40 


75 


67.3 ± 4.7 


8.3 ±9.8 


0.33 


0.9 


38.6 


1.5E-11 


8.9 




0.2-0.3 


1.63 


108 


107.0 ± 5.9 


1.3 ± 12.0 


0.33 


0.1 


32.1 


9.2E-12 


5.6 




0.3-0.4 


2.05 


162 


162.0 ± 7.3 


0.3 ± 14.6 


0.33 


0.0 


37.1 


7.7E-12 


4.7 




0.4-0.5 


2.11 


113 


106.0 ± 5.9 


7.3 ± 12.1 


0.33 


0.6 


42.3 


9.6E-12 


5.8 


SS 


0.5-0.6 


3.40 


231 


222.3 ± 8.5 


8.7 ± 17.5 


0.33 


0.5 


57.1 


8.5E-12 


5.1 




0.6-0.7 


4.93 


370 


361.3± 10.2 


9.1 ±21.8 


0.29 


0.4 


68.3 


5.7E-12 


3.4 




0.7-0.8 


7.13 


649 


600.7 ± 17.1 


49.4 ± 30.7 


0.49 


1.6 


163.9 


8.1E-12 


4.9 




0.8-0.9 


4.69 


331 


310.0 ±8.4 


21.2 ± 20.0 


0.23 


1.1 


86.4 


8.3E-12 


5.1 




0.9-1.0 


2.05 


109 


116.0 ±6.2 


-7.0 ± 12.1 


0.33 


-0.6 


22.1 


5.5E-12 


3.3 




0-1 


30.78 


2216 


2120.3 ±29.5 


96.3 ± 55.5 


0.41 


1.8 


311.6 


4.1E-12 


2.5 




0.0-0.1 


2.89 


511 


509.7 ± 12.9 


1.3 ±26.1 


0.33 


0.0 


102.5 


8.6E-12 


5.2 




0.1-0.2 


3.23 


477 


501.3 ± 12.8 


-24.3 ± 25.3 


0.33 


-0.9 


52.5 


4.6E-12 


2.7 




0.2-0.3 


2.53 


281 


308.0 ± 10.1 


-27.0 ± 19.5 


0.33 


-1.4 


33.6 


4.1E-12 


2.5 




0.3-0.4 


2.53 


266 


230.0 ± 8.7 


36.0 ± 18.5 


0.33 


2.0 


218.1 


2.1E-11 


12.4 




0.4-0.5 


2.79 


248 


261.0 ±9.2 


-13.0 ± 18.3 


0.33 


-0.7 


43.2 


5.4E-12 


3.3 


HS 


0.5-0.6 


2.27 


210 


205.3 ± 8.2 


4.7 ± 16.6 


0.33 


0.3 


73.1 


LIE— 11 


6.8 




0.6-0.7 


2.99 


242 


243.0 ± 8.9 


-1.0 ± 17.9 


0.33 


-0.1 


65.8 


7.1E-12 


4.3 




0.7-0.8 


2.63 


234 


222.7 ± 8.6 


11.3 ± 17.5 


0.33 


0.6 


116.3 


1.1E-11 


6.8 




0.8-0.9 


1.95 


199 


189.7 ±7.9 


9.3 ± 16.2 


0.33 


0.6 


118.5 


1.3E-11 


7.8 




0.9-1.0 


2.13 


235 


261.3 ±9.3 


-26.3 ± 17.9 


0.33 


-1.4 


33.8 


3.9E-12 


2.4 




0-1 


25.9 


2903 


2932.0 ±31.1 


-29.0 ±62.2 


0.33 


-0.5 


191.0 


1.8E-12 


1.1 



Note. — 



Refer to Table[2]for the meaning of the columns. 
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FIG. 8. — Zoom of Figure [2] around the cycle IV campaign between 2008 April 9 and May 2. The open black points in the VHE MAGIC panel show the 
non-significant measured integral fluxes with their statistical error bars (whereas the ULs take into account also the systematic errors). 
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FIG. 9. — Zoom of Figure [2] around the cycle II campaign between 2006 July 14 and August 5. The open black points in the VHE MAGIC panel show the 
non-significant measured integral fluxes with their statistical error bars (whereas the ULs take into account also the systematic errors). 



